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Abstract 

Since a lot of theoretical predictions on the production of new particles (Higgs, 
SUSY) at the Tevatron are based on model estimations of the proton gluon density 
behavior at low x and high values of a transfered momentum Q 2 , the study of a 
possibility of a measurement of the gluon density in this kinematic region directly 
q ■ in Tevatron experiments is obviously of a big interest 0. 

CN . Basing on the selection criteria proposed ealier in [|I|] and |2[], the background 

events suppression factors and corresponding signal events selection efficiencies are 
determined here. The estimation of the number of "7 + jet" events suitable for 
measurement of gluon distribution in different x and Q 2 intervals at Tevatron Run II 
is also done. 

It is shown that with integrated luminosity L int = 3 fb' 1 it would be possible 
1> ', to collect about one million of these events. This number would allow to cover a 
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new kinematical region, 10~ 3 <x< 1.0 with 1.6 • 10 3 < Q 2 < 2 • 10 4 (GeV/c) 2 , not 



O 

O ' studied in any previous experiment. This area includes the values of Q 2 that are, on 

the average, by about one order of magnitude higher than those reached at HERA 
. now. The rates of g c — > 7 dir + jet events are also obtained. 

O I Key-words: gluon distribution, proton, direct photon 

Selection cuts and background suppression. 

To estimate the efficiency of the selection criteria used in 0-0 we carried out the simu- 
lation [] with a mixture of all QCD and SM subprocesses with large cross sections existing 
in PYTHIA [| (namely, in notations of PYTHIA, with ISUB=1, 2, 11-20, 28-31, 53, 
68). The events caused by this set of the subprocesses may give a large background to 
the "7 dir + jet" signal events defined by the subprocesses ([[]) and (fj) []: 
"Compton-like" scattering 

qg->q + i (i) 

and "annihilation" process 

qq-^ g + 1 (2) 

(ISUB=29 and 14) that were also included in the performed simulation. 

Three generations with the above-mentioned set of subprocesses were done. Each of 
them was fulfilled with a different value of the PYTHIA parameter CKIN(3) = p™ m that 
defines the minimal value of parton P t appearing in the final state of a hard 2 — > 2 parton 
level fundamental subprocess in the case of the initial state radiation (ISR) absence. These 
values were p™ 171 = 40, 70 and 100 GeV/c. By 40 million events were generated for each of 

1 PYTHIA 5.7 version with default CTEQ2L parameterization of structure functions is used here. 
2 A contribution of another possible NLO channel gg — > gj (ISUB=115 in PYTHIA) was found to be 
still negligible even at Tevatron energies. 
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the ft™ m value. The cross sections of the above-mentioned subprocesses define the rates of 
the corresponding physical events and, thus, serve during a simulation as weight factors. 

We selected "7 dir -candidate +1 jet" events containing one 7 cfer -candidate (denoted in 
what follows as 7) and one jet (found by a subroutine LUCELL implemented in PYTHIA 
as a jetfinder) with P t 3et > 30 GeV/c. Here and below, as we work at the PYTHIA particle 
level of simulation, speaking about the 7 dir -candidate we actually mean, apart from 7 , 
a set of particles like electrons, bremsstrahlung photons and also photons from neutral 
meson decays that may be registered in one DO calorimeter f] cell of the At] x A<fr = 0.1x0.1 
size. 

Below we consider a set of 17 cuts that are separated into 2 subsets: a set of the 
"photonic" cuts and a set of the "hadronic" ones. The first set consists of 5 cuts used to 
select an isolated photon candidate in some Pt 1 interval. The second one, that includes 12 
cuts, deals mostly with jets and clusters and is used to select events having one "isolated 
jet" and a limited P t activity out of "7 + jet" system. 

All of the cuts are listed in Table [1]. Their influence on the signal-to-background ratio 
Sj B is presented in Table 2 (for a case of the most illustrative intermediate interval of 
event generation with pj" m = 70 GeV/c) and Tables 3-5. 

Tables 1 and 2 are complementary to each other. The numbers in the left-hand column 
( "Cut" ) of Table coincide with the numbers of cuts listed in Table IB. 

The second and third columns of Table 2 contain, respectively, the numbers of signal 
direct photons (S) and background ^ dlT — candidates (B) left in the sample of events after 
application of each cut from Table 1. The numbers of background events B do not include 
events with electrons that fake photon signal. Their numbers in the samples are presented 
separately in the last right-hand column "e ± " . The other columns of Table |2] include the 
values of efficiencies Effg(B) (with their errors) defined as a ratio of a number of the 
signal (background) events that passed under each of the cuts (1-17) to the number of 
the preselected events (1st cut of this table). They are followed by the column containing 
the values of S/ B ratio. 

Table 1: List of the applied cuts (will be used also in Tables Q-§p. 

1. a) > 40 GeV/c, b) P t jet > 30 GeV/c, 9. A</> < 17°; 

c) \rf<\ < 2.5, d) P t hadr < 7 GeV/c*; 10. P t miss /P{<< 0.10; 

2. P t isol < 5 GeV/c, e*t < 15%; 11. P t dust < 20 GeV/c; 

3. Pt" 1 > y™ in \ 12. P t clust < 15 GeV/c; 

4. Ptifnn < 1 GeV/c**; 13. P t clust < 10 GeV/c; 

5. p/ s °?< 2 GeV/c, < 5%; 14. P t out < 20 GeV/c; 

15. P t out < 15 GeV/c; 



6. Njet < 3 

7. Njet < 2 

8. Njet = 1 



16. P t out < 10 GeV/c; 

17. e^ et < 3%. 



* maximal Pt of a hadron in the ECAL cell containing a 7 dir -candidate; 
** A scalar sum of P t in the ring: P t sum (R = 0.4) - P t sum (R = 0.2). 

Line number 1 of Table 1 contains four primary preselection criteria. It includes general 
Pt cuts (la), (lb) as well as the cut connected with the DO electromagnetic calorimeter 
(ECAL) geometry (lc) and the cut (Id) that excludes 7 dir -candidates accompanied by 
hadrons entering the same calorimeter cell as 7. 

Line number 2 of Table 1 fixes the values of 7 dir -candidate isolation parameters P t ls ° l 
and e 7 , i.e. the value of the scalar sum of Pt of all particles surrounding 7 cfer -candidate 



3 The geometry of DO detector was used in the simulation Q. 
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within a cone of R] sol = ((At/) 2 + (A0) 2 ) 1 / 2 = 0.7 and its fraction = P t iso[ / P t \ 

The third cut selects the events with 7 d * r -candidate having P t higher than p™ m thresh- 
old. We impose this cut to select the samples of events with P t 7 > 40, 70 and 100 GeV/c 
as ISR may smear the sharp kinematical cutoff defined by CKIN(3) ||. This cut reflects 
an experimental viewpoint when one is interested in how many events with 7 d * r -candidates 
are contained in some definite interval of Pf. 

The forth cut restricts a value of P™°n g = PtR%A - ^tI=o.2> where P t l ^° l is a sum 
of P t of all ECAL cells contained in the cone of the radius R around the cell fired by 
7 dir -candidate @, @. 

The fifth cut makes tighter the isolation criteria within R = 0.7 than those imposed 
onto 7 cfer -candidate in the second line of Table 1 . 



Table 2: Values of significance and efficiencies for p™ m = 70 GeV/c. 



Cut 


S 


B 


Effs(%) 


Eff B (%) 


S/B 


e± 


1 


39340 


1247005 


100.00± 0.00 


100.000± 0.000 


0.03 


17562 


2 


36611 


51473 


93.06± 0.68 


4.128± 0.019 


0.71 


4402 


3 


29903 


18170 


76.01± 0.58 


1.457± 0.011 


1.65 


2038 


4 


26426 


11458 


67.17± 0.53 


0.919± 0.009 


2.31 


1736 


5 


23830 


7504 


60.57± 0.50 


0.602± 0.007 


3.18 


1568 


6 


23788 


7406 


60.47± 0.50 


0.594± 0.007 


3.21 


1554 


7 


23334 


6780 


59.31± 0.49 


0.544± 0.007 


3.44 


1460 


8 


19386 


4136 


49.28± 0.43 


0.332± 0.005 


4.69 


1142 


9 


18290 


3506 


46.49± 0.42 


0.281± 0.005 


5.22 


796 


10 


18022 


3418 


45.81± 0.41 


0.274± 0.005 


5.27 


210 


11 


15812 


2600 


40.19± 0.38 


0.208± 0.004 


6.08 


176 


12 


13702 


1998 


34.83± 0.35 


0.160± 0.004 


6.86 


130 


13 


10724 


1328 


27.26± 0.30 


0.106± 0.003 


8.08 


88 


14 


10636 


1302 


27.04± 0.30 


0.104± 0.003 


8.17 


86 


15 


10240 


1230 


26.03± 0.29 


0.099± 0.003 


8.33 


84 


16 


8608 


984 


21.88± 0.26 


0.079± 0.003 


8.75 


64 


17 


6266 


622 


15.93± 0.22 


0.050± 0.002 


10.07 


52 



(*) The background (B) does not include the contribution from the "e^ events" (i.e. in which fake 
7-candidate) that is shown separately in the right-hand column "e . 

The cuts considered up to now, apart from general preselection cut P t 3et > 30 GeV/c, 
used in the first line of Table [I], were connected with photon selection ( "photonic" cuts) . 

The cuts 6-9 are connected with the selection of events having only one jet and the 
definition of jet-photon spatial orientation in 0-plane. The 9-th cut selects the events 
with jet and photon transverse momenta being "back-to-back" to each other in 0— plane 
within the angle interval of the A0 = 17° size 

In line 10 we used the cut on p t miss to reduce the background contribution from the 
electroweak subprocesses qg — > q' + W ± and q q' — > g + W ± with the subsequent decay 
W —>■ e v that leads to a substantial p t miss value (see plots in 0]). One can see from 
the last column of Table [2] ll e ±n that the cut on p t miss reduces strongly (in about 4 times) 
the number of events containing as direct photon candidate. 

Moving further we see from Table |2| that the cuts 11-16 of Table [TJ reduce the values 
of transverse momenta of clusters (mini-jets) p t clust and the modulus of a vector sum of 

4 i.e. within the size of three calorimeter cells 
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all particles that are not included into the "7 + jet" system, i.e. P° u , down to the values 
less than 10 GeV/c. The 17-th cut of Table p] imposes the jet isolation requirement. It 
leaves only the events with jets having the sum of P t in the ring of the AR = 0.3 size 
surrounding a jet to be less than 3% of Pt jet . From comparison of the numbers in the 
10-th and 17-th lines we make an important conclusion that all these new cuts (11-17) 
may lead to the following about two-fold improvement of S/ B ratio Q. This improvement 
is achieved due to a reduction of P t activity out of "7 + jet" system |. 

It is also important to mention that the total effect of "hadronic cuts" 6-17 for the case 
of p™ m = 70 GeV I c consists of about twelve-fold decrease of background contribution at 
the cost of less than four-fold loss of signal events (what results in about 3.2 times growth 
of S/B ratio). So, in this sense, we may conclude that from the viewpoint of S/B ratio 
a study of "7 + jet" events may be more preferable as compared with a case of inclusive 
photon production. 



Table 3: Number of signal and background events remained after cuts. 



~ mm 
P± 

(GeV/c) 


Cuts 


7 

direct 


7 

brem 


photons from the mesons 


e± 




V 






40 


Preselected 


18056 


14466 


152927 


56379 


17292 


14318 


2890 


After cuts 


6238 


686 


824 


396 


112 


104 


24 


+ jet isol. 


3094 


264 


338 


150 


40 


44 


14 


70 


Preselected 


39340 


63982 


761926 


269666 


87932 


63499 


17562 


After cuts 


8608 


424 


320 


146 


58 


36 


64 


+ jet isol. 


6266 


262 


206 


90 


40 


24 


52 


100 


Preselected 


56764 


111512 


970710 


346349 


117816 


91416 


38872 


After cuts 


11452 


280 


124 


92 


24 


24 


136 


+ jet isol. 


9672 


204 


92 


64 


24 


20 


120 



Table 4: Efficiency, S/B ratio and significance values in the selected events without jet 
isolation cut. 





s 


B 


Effs(%) 


Eff B (%) 


S/B 


s/Vb 


40 


6238 


2122 


34.55±0.51 


0.831±0.018 


2.9 


135.4 


70 


8608 


984 


21.88 ± 0.26 


0.079 ± 0.003 


8.8 


274.4 


100 


11452 


544 


20.17 ± 0.21 


0.033 ± 0.001 


21.1 


491.0 



Table 5: Efficiency, S/B ratio and significance values in the selected events with jet 
isolation cut. 





S 


B 


Eff S (%) 


Eff B (%) 


S/B 


s/Vb 


40 


3094 


836 


17.14±0.33 


0.327±0.011 


3.7 


107.0 


70 


6266 


622 


15.93 ± 0.22 


0.050 ± 0.002 


10.1 


251.2 


100 


9672 


404 


17.04 ± 0.19 


0.025 ± 0.001 


23.9 


481.2 



Table includes the numbers of signal and background events left in three generated 
event samples after application of cuts 1-16 and 1-17. They are given for all three intervals 
of Pt 1 . The summary of Table ^| is presented in the middle section (^pj am = 70 GeV/c) of 

5 One must keep in mind that the starting value of S/B ratio has a model dependent nature. 
6 See Jl| and Q for a detailed definitions of p t out and jet isolation parameter e- jet . 
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Table |3| where the line "Preselected" corresponds to the cut 1 of Table [1] and, respectively, 
to the line number 1 of Table ^| presented above. The line "After cuts" corresponds to 
the line 16 of Table |2] and line "+jet isolation" corresponds to the line 17 of Table ^. 

Table [3] shows in more detail the origin of 7 <fer -candidates. The numbers in the column 
"7 — direct" correspond to the numbers of signal events caused by QCD subprocesses (1) 
and (2) left in each of Pf intervals after application of the cuts defined in lines 1, 16 and 
17 of Table [l] (see also column "5"' of Table fj). Analogously, the numbers in the column 
"7 — brem" of Table [| correspond to the numbers of events with the photons radiated 
from quarks participating in hard 2^2 parton interactions. Columns 5 - 8 of Table |3| 
illustrate the numbers of the "7 — mes" events with photons originating from 7r°, rj, uj 
and Kg meson decays. In a case of ft™ m = 40 GeV/c the total numbers of background 
events, i.e. a sum over the numbers presented in columns 4 - 8 of Table ^ are shown in 
lines 1, 16 and 17 of column U B" of Table |^. 

The other lines of Table |3| for p™ in = 40 and 100 GeV/c have the meaning analogous 
to that one described above for p ™ m = 70 GeV / c. 

The last column of Tabled shows the number of preselected events with ^ dir — candidates 
faked by e ± . One can see that at pj™ m = 40 GeV jc their contribution at the level of the 
last cut "+jet isol." of Table 3 (or cut 17 of Table 1) to the value of the considered total 
background is about 2% and grows up to 30% at p™ in = 100 GeV jc. Furhter suppression 
of the events with 7 = e^ depends on the value of track finding efficiency. 

Here we take the efficiency of each cut of Table 1 to be equal to 100% as we study the 
results of simulation at the particle level. For estimation of the size of possible detector 
effects one may take into account the value of a track finding efficiency in the central 
region (\rj\ <0.9) of DO detector (§, 0) found in Run I to be about 83%. 

The numbers in Tables |] (without jet isolation cut) and |5] (with jet isolation cut) 
accumulate in a compact form the final information of Tables [X| — So, for example, the 
values in columns S and B for p™ m = 70 GeV j c include the total numbers of the selected 
signal and background events taken correspondingly at the level of 16-th (for Table |]) 
and 17-th (for Table §) cuts from Table Q. 

It is seen from Table |]that the ratio S/B grows from 2.9 to 21.1 while P t J increases 
from Pt< > 40 GeV/c to > 100 GeV/c interval. 

The jet isolation requirement (cut 17 from Table [I]) noticeably improves the situation 
at low Pt 1 (see Table |5|). After application of this criterion the value of S/B increases 
from 2.9 to 3.7 at P^ > 40 GeV/c while at > 100 GeV/c the value of S/B changes 
only from 21.1 to 23.9 due to a tendency of the selected event sample to contain more 
events with an isolated jet as P/ 1 increases (see [|[ and 0). 

Let us underline here that, in contrast to other types of background, the "7 — brem" 
background has an irreducible nature. Thus, the number of "7 — brem" events should be 
carefully estimated for each P/ 1 interval using the particle level of simulation in the frame- 
work of event generators like PYTHIA. The contribution of "7 — brem" background must 
be taken into account in analysis of the experimental data on prompt photon production 
especially for high P t intervals where, as it is seen from Table 3, it becomes a dominant 
one as compared with the contribution from "7 — mes" and u e ±,J events []. 



7 For a careful theoretical study of '7 — brem" background to direct photon production see |n| [Tl| and 
the corresponding references in them. 
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2 Event rates estimation at the Tevatron. 



One of the promising channels for measurement of the proton gluon density, as was shown 

dir 



111 



121, is a high P t direct photon production pp(p) 



7""' + X. The region of high P t , 
reached by UAl @, UA2 [[14]] , CDF |3J and DO |] extends up to P t « 60 GeV/c and 
recently up to P t = 105 GeV/c |§. These data together with the later ones (see references 



in |16|-[|24|]) and recent E706 E3] and UA6 [Z6] results give an opportunity for tuning the 
form of gluon distribution (see fTTH , [pOf| , p7j ). The rates and estimated cross sections of 
inclusive direct photon production at the LHC are given in [12| (see also pBfl). 

Here for the same aim we consider the process pp — > r ) dxr + 1 jet + X defined in the 
leading order by two QCD subprocesses (1) and (2) (for experimental results see 
and references in [[31]]). The analogous study for LHC energy was done in 0] 



and 



The u j dtr + 1 jet" final state is more preferable than inclusive photon production 
^dir _|_ j£ f rom the viewpoint of extraction of information on gluon distribution. Indeed, 
in the case of inclusive direct photon production the cross section is given as an integral 
over partons distribution functions f a {x a ,Q 2 ) (a = quark or gluon), while in the case of 
pp — > r ) d%T + ljet + X for P t jet > 30 GeV/c (i.e. in the region where ll k t " smearing 
effects are not important, see |^TJ) the cross section is expressed directly in terms of these 
distributions (see, for example, p9|): 



da 



dr)idr)2dPt 



where 



Y] x a f a (x a ,Q 2 )x b f b (x b , Q 2 )i- (ab 

a,b dt 



x a ,b = Pt/Vs- (exp(±rji) + exp(±r) 2 )). 



34), 



(3) 



(4) 



We used the following designations above: 771 = rf 1 , r] 2 = if , P t = Pt 7 ; a,b = q, g, g; 3,4 
= q, q, g, 7. Formula (|3|) and the knowledge of the results of independent measurements 
of q, q distributions allow the gluon distribution f g (x, Q 2 ) to be determined after account 
of selection efficiencies of 7 cfar -candidates and the contribution of background, left after 
the used selection cuts (1-13 of Table P, as it was discussed in Section 1 keeping in mind 
this application. 

In the previous works (e.g. see [|lj _ 0) a lot of details connected with the structure and 
topology of these events and the objects appearing in them were discussed. Having all 
this information we are in position to discuss an application of the "7+ jet" event samples 
selected with the proposed cuts and to estimate the rates of gluon-based subprocess (1). 

Table [5] shows the percentage of "Compton-like" subprocess (1) (amounting to 100% 
together with (2)) in the samples of generated with PYTHIA events and then selected 
with cuts 1 - 13 of Table 1 for different Pf and rf^ intervals: Central (CC, |7/| < 0.7), 
Intercryostat (IC, 0.7< |r/| < 1.8) and End (EC, 1.8< \rj\ <2.5) calorimeter. 

Table 6: The percentage of "Compton-like" process q g — » 7 + q. 



Calorimeter 
part 


interval (GeV/c) 


40-50 


50-70 


70-90 


90-140 


CC 


84 


80 


74 


68 


IC 


85 


82 


76 


70 


EC 


89 


85 


82 


73 



In Table [7| we present the Q 2 {= (Pt 1 ) 2 ) and x (defined according to (||)) distribution 
of the number of events (for integrated luminosity Li nt = 3 fb~ l ) that are caused by 
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subprocess (1) and passed cuts 1 - 13 of Table 1 with P t iso1 < 4 GeV/c and e< < 7% 
(applied in the line 5) Q. 

The analogous information for events caused by (1) with the charmed quarks gc —> 
^dir _|_ c j g presented in Table ^. The simulation of the process gb —> r ) dir + b shows 
that the rates for 6-quark are 8-10 times lower than those for c-quark. 



Table 7: Number of g q — > ^ dtr + q events at different Q 2 and x intervals for Li nt = 3 fb 



Q 2 

(GeV/c) 2 


x values of a parton 


All x 


.001 - .005 


.005 - .01 


.01 - .05 


.05- .1 


.1 - .5 


.5-1. 


.001 - I. 


1600-2500 


8582 


56288 


245157 


115870 


203018 


3647 


632563 


2500-4900 


371 


13514 


119305 


64412 


119889 


3196 


320688 


4900-8100 





204 


17865 


13514 


26364 


1059 


59007 


8100-19600 








3838 


5623 


11539 


548 


21549 




1 033 807 



Table 8: Number of g c — > r j dir + c events at different Q 2 and x intervals for Li nt = 3/6 1 . 



Q 2 

(GeV/c) 2 


x values of a parton 


All x 


.001 - .005 


.005 - .01 


.01 - .05 


.05- .1 


.1 - .5 


.5 - 1. 


.001 - 1. 


1600-2500 


264 


2318 


21236 


11758 


14172 


58 


49805 


2500-4900 


13 


332 


9522 


6193 


7785 


40 


23885 


4900-8100 





4 


914 


1055 


1648 


16 


3637 


8100-19600 








142 


329 


612 


8 


1092 




78 419 



Fig. 

covered 
is given 



1 shows in the widely used (x, Q 2 ) kinematic plot (see |32H ) what area can be 
by studying the process q g —* 7 + q. The number of expected events in this area 
in Table [7|. From this figure and Table |7| it becomes clear that at integrated 

luminosity L int = 3 fb' 1 it would be pos- 
sible to study the gluon distribution with 
good statistics of "7 + jet" events in the 
region of small x at Q 2 about one order 
of magnitude higher than now reached at 
HERA 0, [0. It is worth emphasiz- 
ing that extension of the experimentally 
reachable region at the Tevatron to the re- 
gion of lower Q 2 overlapping with the area 
covered by HERA would also be of great 
interest (for instance, it would allow to 
test analytical solutions of DGLAP equa- 
tions that describe Q 2 -evolution of struc- 




ture functions at small values of x \ 35 1 . 



Figure 1: The (x, Q 2 ) kinematic region for pp —>■ r y dir + jet process. 



An application of cuts 14-17, as it is seen from Table 2, leads only to 20% improvement of S/B ratio. 
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